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Abstract 

The finite temperature and finite density dependence of the neutron-proton 
mass difference is analysed in a purely hadronic framework where the p — 
ijj mixing is crucial for this isospin symmetry breakdown. The problem is 
handled within Thermo Field Dynamics. The present results, consistent with 
partial chiral and charge symmetry restoration, improve the experimental 
data fit for the energy difference between mirror nuclei. 
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I. INTRODUCTION 



Symmetry restoration in spontaneously broken gauge theories has been one of the first 
apphcations of field theory at finite temperature [|l|. In the Standard Model, such a phase 
transition is supposed to have taken place in the early Universe when the temperature was 
of the order of a few hundred GeV. Different types of phase transitions are also expected 
to occur at much lower temperatures and with hadronic non zero densities. In connection 
with it, the broken chiral symmetry has been predicted to be restored at about the same 
temperature at which the deconfinement sets in. This temperature is in the range of 150 
to 200 MeV as it is indicated from lattice calculations and phenomenological models 
based on effective Lagrangians High density conditions are supposed to have similar 
consequences on the hadronic matter. 

In the nuclear scale there are also indications of partial deconfinement in order to explain 
the EMC effect 0. It has been argued that nucleons in nuclei occupy a larger volume than 
nucleons in vacuum , a phenomenom that could be understood as a nucleon swelling . It 
has been also analysed in connection with the Nolen Schiffer Anomaly (NSA) [0 for different 
models of the nucleon [||- [|10| . 



In recent years, the possibility of producing quark-gluon plasma by means of very en- 
ergetic nucleus-nucleus collisions []TT| opened a rather important program of investigation 
of matter under extreme conditions that could shed light on the fundamental problems 
mentioned above. 

In the present work we analyse the behavior of the neutron-proton mass difference with 
temperature and density, within meson theory. Our results are consistent with both par- 
tial chiral and charge symmetry restoration, with increasing temperature and density. In 
particular, the density effects are in order to clear away the Nolen-Schiffer Anomaly. 

Although a full description of the neutron-proton mass difference calls for a non pertur- 
bative approach, it is still lacking. Nevertheless, one may gain a good understanding of the 
problem, relying on the perturbative methods used in many body nuclear physics. There, 
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relativistic perturbation theory (RPT) is mainly used for the analysis of hadron interactions 







In a previous article ||13|, we have obtained a good outcome for the n-p mass splitting 
using RPT at the hadronic level. In that work we have shown that the role of the rho-omega 
mixing interaction is crucial in the understanding of the n-p mass difference, in a hadronic 
context (see also |^). We have concluded that the mixing of the vector mesons coming 
from the u-d mass difference, has important consequences on the nucleon self-energy and is 
the main non-electromagnetic charge symmetry breaking (CSB) contribution to be consid- 
ered. This effective interaction has been thoroughly investigated since 1960 fl^, particularly 



within the tadpole picture [0- |T^. It has been generally claimed that the main source for 
that vertex has a non-electromagnetic origin related to the quark mass difference |[T7| , p!8| . 
The connection between the mixing lagrangian of eq.(l) and the microscopic quark level 



has been recently explored within different models: QCD sum rules [19|, Coleman-Glashow 



tadpoles pO|, constituent quark model |18| and the Nambu-Jona-Lasinio model ET^ 



Notice that AM^^L^p = 1.29 MeV while AM^'J^p = -0.66 to -0.76 MeV [0,0, implying 



a 2 MeV strong contribution. Within the scheme of ref. |T3[ this quantity can be obtained 
from the nucleon self-energy given in eq.(2), using the experimental values of couplings and 
meson masses m the literature |5|- [|ll, H. 



II. THE MODEL 

In the present work we obtain the finite temperature and density (FTD) dependence of 
the nucleon mass splitting using the same framework as in ref. [|13[ where we have considered 
the nucleon and the vector mesons as the fundamental dynamical degrees of freedom. In 
this context, the calculation starts from the following Hamiltonian |1^ 



+\g^N{pi)YujM)N{p) + ApJ(g)^'^(g) (1) 
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which corresponds to the standard minimal formulation of the interaction under considera- 
tion. Here gp and g^^ are the vector meson coupling constants, is the isovector anomalous 
magnetic moment, M„ = Mp = M, irtp and are the nucleon and mesons masses respec- 
tively and A is the mixing matrix element. The lowest order self-energy correction to the 
nucleon mass coming from which explicitely breaks the charge symmetry is 

[p) 



(2vr)4 

Y^'Sip - g)|(7, - ^^cr,.qnr3XD4q)Dp{q)\ (2) 



where 



S{p-q) = j^^:^ and iP.(g) = l/(g^ - m^). 

We will also include form factors at the meson-nucleon vertices to take into account the 
effect of the nucleon structure in a phenomenological way, thus giving to the model a wider 



range of application in its momentum dependence ||T3| 



y\2 _ ^2 

Fv{q) = ^ f- where v = p,uj 

Af, - q^ 

A natural framework for the study of matter under FTD conditions is the so-called 
Thermo Field Dynamics (TFD) p5[. TFD is a real time formalism of the statistical field 
theory, very powerful for describing non-isolated many body systems. It is a canonical field 
theory formulation in which the Hilbert space is doubled and each field operator has two 
independent components belonging to the thermal doublet. Correspondingly, the Green's 
functions, self-energies etc., are expressed by the thermal matrices. By virtue of this ex- 
tension of the Hilbert space, the pathologies of the pioneer formulations of the real time 
formalism approach are avoided [^. Moreover because the Gell-Mann-Low formula and 
the Wick's theorem for the perturbation expansion are available in TFD, the usual perturba- 
tion theory at zero temperature and density can be easily extended to FTD. This formalism 



is also particularly useful to perform both high and low temperature expansions, a feature 
much less accesible in the imaginary time formalism Consequently, perturbation theory 
is at hand using the Feynman diagrams technique proper of RPT. Moreover, the TFD free 
propagators can be explicitely separated into two parts: one term being the usual one and 
the other part depending on temperature and density. For the sake of brevity we will not 
discuss the derivation of the corresponding thermal matrices for the modified propagators 
to be used in the present calculation For fermions one gets 



s:w = s:y{k)+ss'{k) 

/ 
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^ sin^Ok,) isin(2^fcj 
^ isin(26'fcj -sin^iOko) 



(3) 



with 



e(fco) 

(l + e--)V2 



e(-fco) 

(l + e-)V2 

e^/^9(-fco) 

(l + e-^)V2 ~ (l + e-)i/2 



(4) 



where z = [ko — fi)/kBT, fi is the relativistic chemical potential and 0(/co) is the step 
function. Similarly for the meson propagator one obtains 



D^\k) = D^^^^i^k) + D^^^'i^k) 

( ^ 1 ^ 

k'^—m'f+ie 







k'^ —rn^—ie j 
( 

- 2ni5{k^ - ml) 



sinh^(0fcj |sinh(20fco) 



V 2 



|sinh(20fcg) — sinh^ 



(5) 



where the trigonometric functions for the particle-antiparticle distributions have to be re- 
placed by the corresponding hyperbohc ones ||2^. The (a, b) = (1, 1) elements of the thermal 
matrices are the standard Feynman physical propagators while the (a, b) = (2, 2) are the 
so-called thermal ghost partenaires and the off-diagonal terms are mixed in nature. TFD 
establishes that inner vertices of a diagram can be of either type, physical or ghost while the 
external ones ought to be physical. In view of this, and up to second order perturbation, only 
the mixing vertex A involve both kind of terms. However, for the ghost type, the product 
of meson propagators will produce a vanishing contribution to the loop. Consequently, the 
full FTD expression for the second order self energy results when S^^ replaces 5* and Dl^ 
replaces in eq.(2). Now we can separate S as 

S = + loop{S''''Df'^Dl}'^) 



where the first term, = loop{S^D^D^), represents the usual zero temperature and 



density Feynman contribution calculated in Ref. [jT3[, and the others are finite T-5 correc- 
tions. On the other hand, for nucleons on-shell 6Mp = Up{p)T,Up{p) = —Un{p)^Un{p) , with 
p"^ = M^, then, performing the go integration, we obtain for the second term in eq.(6) (see 
%(!)) 

Ve^ -M2 ((8 + 6k^)M - (4 + 6k^)e ) 

. 2M2-m2 ^ w 2M2~m2 ^ w 2Af2-A2 ^ w 2Af2-A2 ^ ^ 
V 2M ^ 2M ^ ''v 2M ^ 'y 2M ^ ' 

where the contribution coming from antiparticles has been omitted since we will consider 
temperatures much below pair production, A;_b/?~^ << Mc^. Furthermore, as the presence 
of real mesons calls for extreme conditions, the in-medium corrections to the bosonic prop- 
agators included in the last two terms of eq.(6), shall be neglected. 
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Since our description has been given in terms of effective fields and couplings, it wouldn't 
make sense to explore our system in the high T-5 regime where total deconfinement could 
take place. Hence it is natural to consider T « M and densities not too high. For 
arbitrary T and 5 conditions, eq.(|^) must be solved numerically, nevertheless one can find 
usefuU expressions for some physical scenarios. To this end one has to find accordingly, 
aproximated formulas to the chemical potential which is defined in terms of temperature 
and density by 

^ " {2nn)^ I ^ ^ (e(fco-M)/fcsT e(fco+M)/fcsT + i) 

At finite density, the chemical potential (in units c = h = = 1) can be aproximated by 

/z^M + Tlog^(^P (9) 



provided that [pO[ 

27r(-)2/3 « TM « (10) 

Using eq.(^ at low temperatures gives a very small value for the self-mass correction S'^Mp 
(see next section). In particular, vanishes as T ^ 0. 

On the other hand, in the range of the nuclear matter density, a computation of n 
results in 



which is a good aproximation for 



9^4/3 

T' « TT-rdnr-TTT. (12) 



(37r252/3 + ]Vf2) 

Now, performing a standard low T expansion in eq.(^, we get 
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Mo 



2/i2 - M2 

6^ ^-^-"^ -^-"^(/ig-M2)/.g' 



V-ng\^,,) - g{^o) ,T 2 )) (13) 



with 



-M2 ((8 + Qk^)M - (4 + 6A;^)e 



2M 



)( 



2M2- 



2M 



-e){ 



2M2-A2 
2M 



)( 



2Af2-A 2 
2M 



(14) 
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In this context, it can be seen from eq.(13), that is significant even at T = 0. Equa- 

tions ([Tl| ) - ([T3|) are good aproximations for astrophysical temperatures (a couple of MeV's) 
and baryon densities in the range of the nuclear matter density 5q = 0.1934/m~^ [STj. 
Under these conditions, the T^-dependent term of eq . ([T3|) is negligible. Hence, the main 
contribution comes from the first term, which only depends on 6 (see next section). 

Therefore, at finite density equations (^ and ( pT]) have smooth limits for T — ^ 0, and 
so does eq.(0). On the other hand, at finite temperature and very low densities, one has 
fi = 6 ■ f{j3), implying that the limit 5 — is also smooth for the self-mass correction 
(0). In this case, S^Mp depends on e~*^/^ and is exponentially small for T « M. The 
analysis above shows that in the limit of vanishing chemical potential and zero temperature, 
S'^Mp — i> as expected. 



III. NUMERICAL RESULTS AND FINAL CONCLUSIONS 

In order to obtain numerical predictions from this approach we take experimental val- 
ues for all the physical quantities appearing in our formulas: m^^ = 783 MeV, nip = 770 
MeV, M = 939 MeV, = 3.7. Concerning the coupling constants we work with a set 
9p9lu^ = 0.492 GeV^ chosen in order to saturate the 2 MeV hadronic contribution to the 
zero temperature and density nucleon mass difference within the model. This election is 
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well inside the accepted range of variation of these coupling constants found in the litera- 



ture [|T5|,|T6|,|T8|,|4[. 

As it is known, eq.(|^) is mostly suitable for high temperatures and low densities, nev- 
ertheless as our model is not expected to be reliable far from deconfinement (To ~ 180 
MeV), we should consider MT below T|, instead of M^, together with conditions (JT^). For 
example, at T < 15 MeV and densities satisfying these conditions (say below (5o/10), the 
value of A^M = (5^M„ — 5^Mp is below a couple of keV, which is negligible (recall that in 
the chosen units, 1/m = 197 MeV^^). For low temperatures and baryon densities ranging 
from 5o/2 to 1.55o, we have to use eqs.(|TTD-(p!3D. In this case, in turn, A'^M goes from —0.1 
to -0.25 MeV (see fig. 2). 

Evidently, it is the non-trivial shape of /i(5, T) which is responsible for the behavior of 
the in-medium mass difference. As we have mentioned, the temperature contribution is 
negligible for astrophysical temperatures, which we are interested in. 

In connection with these results on the neutron-proton mass difference as a function of 
the nuclear density, let us finally comment on the Nolen-Schiffer Anomaly. The NSA is a 
persistent problem in nuclear physics and the anomaly is the failure of theory to explain the 
mass differences between mirror nuclei {i.e nuclei with Z = A ±1/2 and N = 1/2); 
a gap amounting to a few hundred keV. The effects of nuclear structure have been widely 
discussed with only partial success . The mass difference of mirror nuclei can be written 
as Mz> — Mz< = '^Eem — {Mn — Mp) where M^> is the mass of the nucleus with the larger 
charge, AEem is the electromagnetic self energy difference between the nuclei and M„ — Mp 
is the nucleon mass difference inside the nucleus. 

Since AEem has been exhaustively analysed, in recent years particular attention has been 
paid to the second term. A variety of models have been put forward in order to avoid this 
problem. Generally it is found that M„ — Mp is a decreasing function of the nuclear density 
0- [|T^. It means that within these models high density is also expected to produce a partial 
charge symmetry restoration. This is the expected behavior to deal with the anomaly. 

In this way, the final results of our calculation are in the right direction to remove the 

8 



anomaly. 

It should be mentioned that some recent literature has suggested that the A (off-shell) 



momentum dependence could be important |p5[. However, this conclusion resulted from 



the use of oversimplyfied models for the p — uj correlation function. In fact, it has been 
recently shown |TP[ that the effect of the different widths of the p and uj mesons implies 
that A is almost momentum independent. This is consistent with the present understanding 
of the CSB phenomenology [^,^], where the p — uj mixing has shown to be of particular 
importance. 

Following similar steps to include the electromagnetic interaction in the current scheme, 
we have found that the nuclear medium has a relatively strong effect on the electromagnetic 
self-energy AMj_p in the opposite direction, ammounting to about 15% of its effect on 
AM^-p- This is in good agreement with the results of Ref. |^ which have been derived 
making a quite different treatment of the problem, in order to include the external conditions 
by means of Skyrme type models. 

In table 1 we show the values of NSA reported in the literature. /^From this table it 
emerges (a rough) accordance between the predictions of the different models (with a not 
very clear density dependent NSA). For these nuclei {6av. — 5o/2) our estimate of the NSA is 
of 0.10 MeV. This is consistent with the values quoted in table 1 both in sign and magnitude 
although not big enough to completely remove the anomaly. 

These results for the n-p mass difference within dense matter, suggest that in this case 
one should include, besides the p — uj mixing, other contributions which, although of minor 
importance in vacuum, seem to be relevant inside nuclei. 



According to the vigor recently regained by the p-u mixing effective interaction [jT9 
in this work we have extended our previous analysis on the mass splitting of an isolated 
nucleon |T^. In the present article we have gone a step further, considering the effects of 
temperature and density on this CSB outcome. 

For standard astrophysical temperatures, far from the deconfining phase, we have found 
a negligible temperature effect on the nucleon self-masses. On the other hand, the density 



effects are significant and were sfiown to produce tfie expected trends to remove tfie NSA. 
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TABLES 



Table I. Values of the anomaly NSA in MeV reported in Refs. []7|, |32| , |3^ ,|9[] . A is the mass 
number of mirror nuclei. 



A 


Nolen-Schiffer 

NSA 


Shlomo 

NSA |32| 


Sato 

NSA [34| 


Hatsuda et al 

NSA g 


0-N 15 




0.16 


0.29 


0.53 


F-0 17 


0.31 


0.31 


0.19 


0.29 


Si-P 29 


0.20 




0.24 




S-Cl 33 


0.24 




0.28 




Ca-K 39 




0.22 


0.43 


0.57 


Sc-Ca 41 


0.62 


0.59 


0.35 


0.42 
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FIGURES 



Fig. 1. Diagrams contributing to the second term in eq. 6. The broken hne represents the 
T — 5 correction to the nucleon propagator, i.e : S^^'^ (see eq. 3). 

Fig. 2. Neutron-proton mass difference as a function of the hadronic density S in units of 
Sq, the nuclear matter density, at T=0. 
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